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Evolution of Population I stars with initial masses ranging within 7OM < Mzams — 13OM is 
considered under various assumptions on the mass loss rate M. The mass-luminosity relation of 
W-R stars is shown to be most sensitive to the mass loss rate M^ a during the helium burning 
phase. Together with the mass-luminosity relation obtained for all evolutionary sequences 
several more exact relations are determined for the constant ratios 0.5 < f% a < 3, where f% a = 
M/M^a- Evolutionary models of W-R stars were used as initial conditions in hydrodynamic 
computations of radial nonlinear stellar oscillations. The oscillation amplitude is larger in W-R 
stars with smaller initial mass Mzams or with lower mass loss rate M due to higher surface 
abundances of carbon and oxygen. In the evolving W-R star the oscilation amplitude decreases 
with decreasing stellar mass M and for M < 1OM the sufficiently small nonlinear effects allow 
us to calculate the integral of the mechanical work W done over the pulsation cycle in each 
mass zone of the hydrodynamical model. The only positive maximum on the radial dependence 
of W is in the layers with temperature of T ~ 2 x 10 5 K where oscillations are excited by the 
iron Z-bump K-mechanism. Radial oscillations of W-R stars with mass of M > 1OM are 
shown to be also excited by the K-mechanism but the instability driving zone is at the bottom 
of the envelope and pulsation motions exist in the form of nonlinear running waves propagating 
outward from the inner layers of the envelope. 
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Introduction 

During the core helium burning the structure of Population I stars with initial mass of Mzams > 
3OM weakly depends on the radial distribution of the mean molecular weight fi. This is due 
to the fact that radiation dominates in the internal energy, so that the hydrostatic equilibrium 
is governed by the gravity and by the gradient of radiation pressure whereas the gas pressure 
is negligible. Moreover, the opacity is mostly due to the Tomphson scattering and radiative 
equilibrium above the convective core is almost independent of the radial distribution of /j,. 
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Both these effects are strongest in stars with high effective temperature (T cff > 5 x 10 4 K) and 
lead to the mass-luminosity relation of W-R stars (Langer, 1989; M. Beech, R. Mitalas, 1992). 

Together with the mass-luminosity relation (Maeder, 1983; Doom et al., 1986; Maeder, 
1987, Maeder, Meynet, 1987) W-R stars obey the mass-radius relation (Langer, 1989; Schaerer, 
Maeder, 1992) and therefore the sound travel time between the stellar center and the stellar 
surface can be considered approximately as a function of the stellar mass M. In recent years 
W-R stars were found to be unstable against radial oscillations (Glatzel et al., 1999; Fadeyev, 
Novikova, 2003; 2004; Dorfi et al., 2006; Fadeyev, 2007; Fokin, Tutukov, 2007; Fadeyev, 2008), 
so that one might assume that evolution of W-R stars is accompanined by decrease of the 
period of their radial oscillations II. However one should bear in mind that if the pulsational 
instability is due to the iron Z-bump /t-mechanism (Dorfi et al., 2006; Fadeyev, 2008) then 
both the amplitude growth rate and the amplitude of radial oscillations depend on the chemical 
composition of the outer layers of the star. Such a dependence is still unclear and in the present 
paper we consider the pulsational properties of W-R stars with various surface abundances of 
helium, carbon and oxygen. 

The chemical composition of outer layers of the evolving star depends on the both initial 
mass Mzams and mass loss rate M. In this paper we consider the Population I stars with initial 
masses ranging within 70 M < Mzams < 130M Q . To evaluate the effect of uncertainties in 
empirical estimates of mass loss rates the evolutionary computations were carried out under 
various assumptions on M. Some stellar models corresponding to the W-R phase were used 
as initial conditions in hydrodynamic computations of nonlinear radial stellar oscillations. In 
the present paper we continue our previous studies (Fadeyev, 2007; 2008) but now we use 
more recent data on the mass loss rate during helium burning that are most important for the 
structure of W-R stars. 

Method of computations 

Calculations of stellar evolution and stellar pulsations were done with methods that in general 
were described in the previous paper (Fadeyev, 2007), so that below only a few important 
changes are noted. 

The mass loss rate during hydrogen burning and at the initial phase of helium burning 
Mh was calculated according to Nieuwenhuijzen and de Jager (1990). The mass loss rate at 
the later stages of evolution M 3q was calculated using the empirical formula by Nugis and 
Lamers (2000) which takes into account both clumping of the stellar wind and dependence of 
the mass loss rate on surface abundances of helium and heavier elements. In computations of 
stellar evolutuion the transition from Mh to M 3a was done during the initial stage of helium 
burning when the effective temperature of the contracting star reaches T eff = 10 4 K. At this 



point of the evolutionary track the difference between M H and M 3a is about a few percent. 
For example, for M Z ams = 8OM and M = M H the star reaches the effective temperature 
T e g = 10 4 K when its mass is about M = 37M and surface abundances of hydrogen and 
helium are X( 1 H) rs 0.1 and X( 4 He) ~ 0.88, respectively. The mass loss rates evaluated 
with formulae by Nieuwenhuijzen and de Jager (1990) and by Nugis and Lamers (2000) are 
Mh = 9.3 x lO~ 5 M /yr and M 3a = 9.6 x lO~ 5 M /yr. It should be noted that during such a 
short phase of contraction of the star the evolutionary track crosses the H-R diagram so fast 
that variation of the threshold effective temperature by a factor of two around T eff = 10 4 K 
does not affect perceptibly the later evolution of the stars. 

Thus, in calculations of stellar evolution the mass loss rate was determined as 



M = I 



/hMh, (hydrogen burning and initial helium burning), 
f3aM 3a , (helium burning at T eff > 10 4 K), 



where constant factors 1 < /h < 2 0.5 < fs a < 3 were introduced in order to evaluate the 
dependence of the results of evolution calculations on uncertainties in mass loss rates. 

In previous papers (Fadeyev, 2007; 2008) the thermodynamic functions of the gas with 
temperature of T < 10 7 K were calculated using the OPAL equation of state data (Rogers 
et al., 1996). Unfortunately, these data are not quite appropriate at late stages of evolution 
when abundances of elements heavier than helium substantially increase in the outer layers 
of the star. In the present study we computed the tables of thermodynamic quantities for 
about three dozen compositions, so that both stellar evolution and pulsational instability can 
be selfconsistently calculated up to the core helium exhaustion when the stellar matter consists 
mostly of carbon and oxygen. Computing the tables of thermodynamic quantities we assumed 
that heavy elements are carbon, nytrogen, oxygen and neon. ZAMS abundances of these 
elements were taken from Rodgers et al. (1996). 

The thermonuclear reaction network was extended to several dozen isotopes from hydrogen 
1 H to nickel 56 Ni. However in comparison with previous computations the larger number of 
reactions did not affect significantly the results of evolution calculations. This is due to the 
fact that evolutionary calculations were terminated just after the helium exhaustion when the 
energy generation is due mostly to reactions of the tripple-alpha process and a( 12 C, 7) 16 0. 



Mass— luminosity relation 

Domination of radiation in the helium burning phase leads to convergence of the evolutionary 
tracks of massive stars on the H-R diagram. In the upper panel of Fig. 1 are shown the parts 
of two evolutionary tracks with initial masses M Z ams = 8OM and M Z ams = 13OM , both of 
them being computed with / H = = 1. At luminosity logL/L = 5.9 masses and effective 




Figure 1: Upper panel: the parts of evolutionary tracks on the H-R diagram for initial masses 
M Z ams = 13OM (solid line) and M ZA ms = 8OM (dashed line) with f u = f 3a = 1. The circles 
indicate the models with central helium abundance Y c = 0.9 and Y c = 0.1. Lower panel: The parts 
of evolutionary tracks of the star Mzams = 1OOM with /h = fz a = 1 (solid line), /h = 2, f 3a = 1 
(dashed line) and /h = 1, /3a = 2 (dotted line). 

temperatures of these stars (M xs 25M , T ef[ ps 8.2 x 10 4 K) differ from one another by 2% and 
0.1%, respectively. 

W-R stars with close values of the stellar mass M, luminosity L and radius R but with 
different initial masses Mzams have significantly different radial distributions of the mean molec- 
ular weight. This is illustrated in the upper panel of Fig. 2 where for two models of W-R stars 
with mass of M — 25M (M Z ams — 8OM and M Z ams — 13OM ) the radial dependencies of 
the helium abundance X( 4 He) are plotted as a functions of the Lagrangean mass coordinate 
M r . The lower helium abundance and the higher carbon abundance in outer layers of the star 
with larger initial mass Mzams is due to the more efficient helium burning during the preceding 
stellar evolution. Homogeneity of the chemical composition within the large mass fraction of 
the star is due to the large extention of convective cores. 

Dependence of the evolution of W-R stars on the mass loss rate is illustrated in the lower 
panel of Fig. 1 where three evolutionary tracks of the star with initial mass of M Z ams = 1OOM 
are plotted for three different cases of the mass loss: (/h = fz a = 1), (/h = 2, f 3a = 1) and 
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Figure 2: Relative mass abundance of helium X( 4 He) as a function of the Lagrangean mass coordinate 
M r in W-R stars with mass of M = 25M . Upper panel: the case of mass loss /h = fs a = 1- Numbers 
at the curves indicate the initial stellar mass Mzams- Lower panel: three different cases of mass loss 
for Mzams = 1OOM . 

(/h — 1) /3a — 2). Doubling the mass loss rate Mh during hydrogen burning leads to decrease 
of the luminosity of the W-R star with mass of M — 25M© by about ~ 3%, whereas the core 
helium abundance increases by AY C ps 0.1. 

The most important role in the structure of W-R stars belongs to the mass loss rate during 
the helium burning. For example, doubling M 3a leads to decrease of the luminosity in the W-R 
star with mass of M — 25M by 9% and at bthe same time to increase of the core helium 
abundance by AY C ps 0.25. Dependence of the radial distribution of helium abundance X( 4 He) 
on the mass loss rate M 3a is shown in the lower panel of Fig. 2. Moreover, increase of M 3a is 
accompanied by the slower growth of the central temperature, so that the end of core helium 
exhaustion occurs at the lower stellar mass and on the H-R diagram the evolutionary track 
extends to lower luminosities. 

Convergence of the evolutionary tracks on the H-R diagram implies the existence of the 
correlation between the stellar mass M and the stellar luminosity L. Fig. 3 shows the mass- 
luminosity diagram with two pairs of evolutionary tracks for two initial stellar masses (M Z ams — 
9OM , M Z ams = 13OM ) and two cases of mass loss: (/ H = 1, f 3a = 1;/h = 2, f 3a = 1). At 
the final point of each track plotted in Fig. 3 the central helium abundance is Y c ps 10~ 3 . It is 
clearly seen that during the significant part of the helium burning phase the stellar mass and 
luminosity are nearly related by the power dependence. The position of the initial point of the 
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Figure 3: Mass-luminosity relation for /h = 1 and /h = 2 at initial masses Mzams = 13OM (solid 
and dashed lines) and Mzams = 9OM (dotted and dashed-dotted lines). For all models fs a = 1. 

power dependence changes with initial stellar mass and mass loss rate. For example, the central 
helium abundance at the initial point ranges from Y c = 0.85 for Mzams — 7OM to Y c = 0.94 
for Mzams — 13OM . For smaller M Z ams or larger M H the helium burning occurs at the lower 
central temperature, so that on the mass-luminosity diagram the part of the track with power 
dependence moves to smaller values of M and L. 

The parts of evolutionary tracks that can be approximated by the power dependence on 
the mass-luminosity diagram are displayed in Fig. 4. The maximum deviation of each track 
from the power dependence is Al = max|AlogL| w 0.01 for 0.5 < f^, a < 2 and Al ~ 0.03 for 
/3a = 3. The linear fit on the (log M, log L) plane for models with 7OM < Mzams < 13OM , 
1 < /h < 2 and 0.5 < f^ a < 3 shown in Fig. 4 is given by 

log L/L Q = 3.675 + 1.568 log M/M . (2) 

The maximum deviation of the parts of evolutionary tracks from relation (j2J) is Al < 0.11 
and its large value is due to the dependence of the coefficients of this relation on the mass loss 
rate M^ a . Therefore, more exact approximation can be obtained for fixed values of f^ a : 

A L < 0.06), 
A L < 0.03), 

(3) 

A L < 0.02), 
A L < 0.04). 

It should be noted that the mass interval within which formulae ([3]) can be applied depends 
on the initial mass Mzams as well as on parameters /h and f^ a . The stellar mass M and the 
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Figure 4: Mass-luminosity relation of W-R stars with 7OM < M ZA ms < 13OM , 1 < / H < 2 and 
0.5 < fsa < 3 The dashed line shows relation (refwrmll). 

central helium abundance Y c corresponding to limits of these intervals are given in Table 1. 
The last column of Table 1 gives the evolution time t ev for this interval. 

When one considers the stellar mass M as a parameter one should bear in mind that 
the stellar evolution becomes slower as M decreases. In Fig. 5 for M Z ams — 1OOM and 
/h = 1 are shown the plots of the life time At ev within the mass interval [M — AM, M], where 
AM = 0.5M Q . The plots shown in Fig. 5 weakly depend on Mzams and /h- The growth of 
the life time At cv with increasing f^ a is due to the slower increase of the central temperature 
and slower conversion of helium into carbon. For initial stellar mass Mzams — 100M Q and 
mass loss rate during hydrogen burning / H = 1 the second half of the life time correspond to 
W-R masses M < 33M , M < 25M , M < 12M Q and M < 7M Q for f 3a = 0.5, 1, 2 and 3, 
respectively. 

Nonlinear radial oscillations 

During evolution of the W-R star the relative radius of the stellar core increases, so that 
both the sound travel time from the center to the surface and the period of radial oscillations 
gradually decrease. At the same time the instability excitation zone moves closer to the stellar 
surface and the amplitude of oscillations decreases due to diminishing mass of the driving zone. 
Decrease of the amplitude and the period of radial oscillations in evolving W-R stars is shown 
in Fig. 6 where for stars with initial masses Mzams — 9OM and Mzams — 12OM are given 
the plots of the maximum expansion velocity of the outer boundary of the hydrodynamic model 
U max expressed in units of the local escape velocity t> csc and the period of radial oscillations IT. 
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relations ([3]). 
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Figure 5: The life time At cv of W-R stars within the mass interval AM = O.5M as a function of 
the stellar mass M. The initial mass is Mzams = WOMq and /h = 1. The numbers at the curves 
indicate correspondig values of f^ a . 




Figure 6: The ratio of the maximum expansion velocity of the outer boundary U max to the local 
escape velocity -u esc (upper panel) and the period of radial oscillations IT in days (lower panel) as a 
function of the stellar mass M. In solid and dashed lines are shown the evoltuionary sequences with 
Mzams = 12OM and Mzams = 9OM , respectively. 

The growth time of radial oscillations is of the order of the stellar dynamic time scale, so 
that oscillations of W-R stars are stronly nonadiabatic. In contrast to many other radially 
pulsating stars oscillations of W-R stars cannot be described in terms of the standing wave 
since the kinetic energy of the pulsating envelope only once per period passes the minimum 
and the maximum. Radial oscillations of W-R stars should be considered as nonlinear running 
waves propagating from the envelope bottom to the outer boundary. That is why in W-R 
stars with mass of M > 15M Q the pulsation constant (Q > 0.1 day) is significantly larger in 
comparison with pulsation constants of stars radially pulsating in the form of standing waves. 
The only exception is W-R stars with mass of M < 1OM where the small-amplitude radial 
oscillations can be approximately represented by nonadiabatic standing waves. 

In the stellar mass range of 4.5M < M < 2OM the pulsation constant of W-R stars can 
be approximately written as 

logQ = -2.6 + 0.1M©. (4) 

Properties of some hydrodynamic models are listed in Table 2. The mass of outer pulsating 
layers is negligible in comparison with the total mass of the star, so that the abundances of 
helium X( 4 He), carbon X( 12 C) and oxygen X( 16 0) are constant through the envelope. The 
mean pulsation period II was evaluated using the discrete Fourier transform of the kinetic 
energy of the oscillating envelope within the time interval 10 2 < t/U < 10 3 . However, strictly 
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Figure 7: The power spectrum of velocity U (u>) in the outer (upper panel) and the inner (lower panel) 
layers of the W-R star with mass M = 16M . 

speaking, the definition of the period of radial oscillations can be applied to W-R stars with 
rather low stellar masses because for M > 15M the frequency of oscillations at the bottom of 
the envelope becomes somewhat higher than that of the outer layers. In Fig. 7 are shown the 
power spectra of the velocity of gas in the outer (r pa R) and the inner (r pa 0.79R) layers of the 
envelope of the W-R star with mass of M — 16M . Contribution of short-period oscillations in 
the inner layers becomes significant in stars with M > 2OM because they affect the radiative 
flux emerging from the outer boundary and the period of light variatins becomes shorter than 
that of hydrodynamic motions in the outer layers of the pulsating envelope. 

As is seen in Table 2 the pulsational properties of W-R stars depend mostly on the mass 
loss rate M 3q whereas the role of the both initial stellar mass M ZA ms and mass loss rate during 
hydrogen burning M H is significantly weaker. This is due to the fact that variations of M 3a are 
accompanied by significant changes of the both stellar luminosity and surface abundances of 
carbon and oxygen. In particular, increase of the mass loss rate during helum burning leads to 
the smaller stellar luminosity and therefore to the smaller nonadiabaticity of stellar pulsations. 

An important consequence of nonlinear radial stellar oscillations is the increase of the mean 
radius of pulsating layers of the gas. This effect is illustrated by Table 2 where the last column 
gives the mean radius of the outer boundary of the hydrodynamic model R expressed in units 




Table 2: Properties of some hydrodynamic models of W-R stars. 
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of the initial equilibrium radius R. 



Mechanism of pulsational instability 

The elementary spherical layer of gas contributing into excitation of pulsational instability per- 
forms the positive mechanical work during the pulsation cycle, that is the integral of mechanical 
work is positive : W — <j> PdV > 0, where P is the total pressure and V is the specific volume. 
Unfortunately, for hydrodynamical models of W-R stars exact calculation of the radial depen- 
dence of the mechanic work W is impossible because of strongly nonlinear radial oscillations. 
The only exception is W-R stars with mass of M < 1OM where pulsation motions are charac- 
terized by a good repetition of pulsation cycles. Formation of W-R stars with so small masses 
implies rather high mass loss rates during the core helium burning phase (fs a > 3). Omitting 
discussion on possibility of the existence of such W-R stars we consider their pulsational prop- 
erties because the results obtained can be generalized to more massive W-R stars undergoing 
substantially smaller mass loss. 

The upper panel of Fig. 8 shows radial dependences of the inetgral of mecahnical work W 
for two models of W-R stars with mass of M — 6M and M = 8M . It is clearly seen that 
excitation of oscillations (W > 0) occurs in the outer layers of the stars and the maximum of 
W moves to the surface as the stellar mass M decreases. To clarify the origin of the pulsational 
instability one should compare the radial dependence of W with that of the amplitude of 
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Figure 8: Radial dependences of the mechanical work over the pulsation cycle W (upper panel) and 
the normalized spectral density of luminosity L r {u) (lower panel) for models of W-R stars with masses 
M = (i.U. and M = 8M . 

radiative flux variations. To this end we consider the spectral density of luminosity 



at the angular frequency uj = 2n/H. The spectral density L r (u) was computed in all La- 
grangeam mass zones using the discrete Fourier transform within time intervals t/U. < 10 3 . 
The lower panel of Fig. 8 shows two radial dependences of L r (u) that are normalized to the 
surface value. 

Coincidence of radial coordinates of maxima of W and L t (uj) (see Fig. 8) allows us to 
conclude that excitation of oscillations is due to the interaction of radiative flux with gas of the 
envelope. In particular, for W > it is necessary that the gas absorbs radiation at maximum 
compression and becomes more transparent at maximum expansion. Fig. 9 shows the plots 
of variations of the gas density p and opacity k in the mass zone with maximum mechanical 
work W of the model of W-R stars with mass M = 8M Q . It is clearly seen that the positive 
mechanical work W is due to the /t-mechanism. The average temperature of gas in this zone 
is T ~ 2 x 10 5 K and the positive temperature derivative of opacity (d\nK/d\nT) p > is due 
to the iron Z-bump. 

Fig. 10 shows the plots of radial dependencies of the gas density p and the spectral density 
of luminosity L r (uo) for outer layers of W-R stars with masses 12M < M < 18M . The 
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Figure 9: Variations of the gas density p and opacity k in the vicinity of the maximum of mechanical 
work W in the model of the W-R star with mass M = 8M . The plots of p and k are arbitrarily 
shifted along the vertical axis. 

narrow maximum of L r (u) (see the lower panel of Fig. 10) corresponds to layers with agerage 
temperature T ~ 2 x 10 5 K where the gas density p and opacity k reach the maximum simul- 
taneously, that is excitation of oscillations is also due to the /t-mechanism. In layers above the 
excitation zone dL r (uj) / dr m 0, that is there are neither excitation nor damping zones. More 
massive W-R stars possess more extended envelopes with smaller gradient of the gas density. 
The excitation zone is at the bottom of the envelope and radial oscillations exist as successive 
nonlinear waves propagating from the outer boundary of the core to the stellar surface. 

Conclusion 

As follows from results presented above the mass and luminosity of Population I massive stars 
are related by the power dependence during the large part of the helium burning phase. Bound- 
aries of the segment of the evolutionary track within which the coefficients of the power de- 
pendence are constant depends on the initial stellar mass Mzams and the mass loss rate M 
during the preceding evolution. At the initial point of the power dependence the central helium 
abundance ranges from Y c « 0.78 (M ZAM s = 7OM , f 3a = 0.5) to Y c « 0.96 (M ZAM s = 13OM , 
fza = 3), whereas the effective temperature is in the range from T c g rs 3 x 10 4 K (Mzams = 
12OM ) to T eff « 5 x 10 4 K (M ZAM s = 70 M ). At the end of the power dependence the central 
helium abundance ranges within 0.01 < Y c < 0.15. Thus, the limits of applicability of the 
mass-luminosity relation depend on the both initial mass Mzams an d mass loss rate M, so that 
relations (j3J) should be used with corresponding boundary values of the stellar mass M listed 
in Table 1. 
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Figure 10: Radial dependencies of gas density p (upper panel) and normalized spectral density of 
luminosity L r (w) (lower panel) in W-R stars with mass M = 12M , 15M and 18M Q . 

In the present study we considered pulsational properties of W-R stars as a function of 
three parameters: the initial stellar mass M Z ams> the mass loss rate during hydrogen burning 
Mh and the mass loss rate during the helium burning phase M 3a . However the most important 
is M^a since its variations lead to significant changes in surface abundances of helium, carbon 
and oxygen. In particular, at lower mass loss rate M 3a radial oscillations of W-R stars have 
larger amplitudes due to higher surface abundances of carbon and oxygen. 

As is seen from Fig. 6 and Table 2 the period of radial oscillations is a sensitive indicator of 
the stellar mass and during evoulution of the W-R star decreases from as 5.5 hr at M = 22M Q 
to ~ 2.6 min at M — 5M Q . Thus, observational estimates of pulsation periods could provide 
with a direct evaluation of the masses of W-R stars. For example, aAccording to Veen et al. 
(2002a, 2002b, 2002c) the pulsation period of WR46 is II « 0.14 day and as follows from our 
hydrodynamical calculations the stellar mass is M xs 20 M Q . 

The present study is confined to hydrodynamical models of W-R stars near their final stage 
of evolution and such a choice are is to their longer life time. However of great interest are 
more massive W-R stars with higher luminosity and much stronger instability similar to that 
of LBV stars. The study of these objects will be presented in the forthcoming paper. 
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